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Abstract 



Universality properties of the distribution of the generalized eigen- 
values of a pencil of random Hankel matrices, arising in the so- 
lution of the exponential interpolation problem of a complex dis- 
crete stationary process, are proved under the assumption that 
every finite set of random variables of the process have a multi- 
variate spherical distribution. An integral representation of the 
condensed density of the generalized eigenvalues is also derived. 
The asymptotic behavior of this function turns out to depend 
only on stationarity and not on the specific distribution of the 
process. 

Key words: complex moments; Pade approximants; random poly- 
nomials. 

MSC 2000: 15B52, 60B20, 62E15 



1 



Introduction 



Let us consider the following moment problem: to compute the 
complex measure defined on a compact set D C W by 

S(*) = E 7,^-0), Gemt(D), Cj^ChVj^h, ^e(T 

.7 = 1 

from its complex moments 

Sk = y z k S(z)dz = y (x+iyY S(x J t-iy)dxdy, k = 0, . . . , n— 1, n = 2p 

D D 

It turns out that 



(1) 



.7 = 1 



The problem is thus equivalent to solve the complex exponential 
interpolation problem for the data s^, k = 0, . . . , n — 1. It is well 
known that the Q, j = 1 , . . . , p are the generalized eigenvalues 
of the Hankel pencil P = [Ui(s), Uq(s)] where 

U (s) = U (s , . . . , s n _ 2 ), Ui{s) = U (si, . . . , s„_i) 



and 



U(s , . . . , s n _ 2 ) = 



S Si . . . Sp-i 



si s 2 



Sp—1 Sp 



S n -2 



Moreover 7j are related to the generalized eigenvector Uj of P by 
7.7 = wJ[s ,...,s p _i] T . 

Equivalently, (j, j = 1, • • • ,P are the roots of the polynomial in 
the variable 2 

cte*[£/i(s) - 2t/o(s)] 
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which is the denominator of the Pade approximant [p — l,p](z) 
to the Z-transform of s k defined by 

oo 

f(z) = £ s k z-\ 

fc=0 



Denoting random variables by bold characters, let us assume now 
that we know an even number n > 2p of noisy complex moments 

8Lk = s k + v k , k = 0, 1,2, ...,n- 1 

where v k is a stationary discrete process, and we want to esti- 
mate S(z) - i.e. p, (7j, Q, j = 1, . . . ,p - from {a fc } fe=0 ,...,„-i. This 
is a well known difficult ill-posed problem which is central in 
many disciplines and appears in the literature in different forms 
and contexts (see e.g. |21171l9 1fTllfT51ITB] ) . All the quantities defined 



above become random. It is therefore relevant to study the distri- 
bution of these quantities under suitable hypotheses on the noise 
affecting the data. 

In this work the case s k = 0, \/k will be considered assuming that 
the noise is represented by a discrete stationary process, white 
or colored. For example the distribution of every finite set of r.v. 
of the process can be multivariate a-stable which is a class of 
distributions closed with respect to addition (up to scale and lo- 
cation parameters), a property consistent with the naive concept 
of noise. Most of the properties proved in the following will require 
that the distribution of every finite set of r.v. of the process is 
spherical and the density function exists ( fT3l Sec. 1.5], [8, Sec. 2]). 
We remember that under this assumption white noise is neces- 
sarily Gaussian. Therefore e.g. a-stable, non-Gaussian (a / 2), 
spherical processes are colored. 

Motivations to consider the pure noise case are twofold. The prob- 
lem of identifying the presence of a signal in a large noise envi- 
ronment arises in many applied contexts. The distribution of the 
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generalized eigenvalues Cj of the random pencil 

P = [U(a h . . . ,a n „i), U(a , . . . ,a n _ 2 )] 



in the pure noise case is a reference to detect the presence of a 
signal. More generally, when solving the noise filtering problem, 
information on the generalized spectral properties of the noise 
are required. In the specific class of complex exponential models 
a generalized spectrum is given by the condensed density ([1]) of 
the generalized eigenvalues ^ of the random pencil P or, equiva- 
lently, of the poles of the [p— l,p](z) random Pade approximant 
to the Z-transform (formal random series) 



In pi] it was proved that if {&k} is a complex Gaussian white noise, 
expressing the condensed density in polar coordinates, the radial 
density weakly converges to a Dirac distribution centered in 1 and 
the phase density is uniform in (0,27r]. Moreover it was proved 
that for n = 2 the condensed density is the uniform measure on 
the Riemann sphere 



In [5] it was conjectured, on the basis of numerical experiments, 
that the condensed density has an universal behavior i.e., in polar 
coordinates, the radial density is Lorentzian, centered in 1 with 
width dependent on n and the phase density is uniform in (0, 2tt}. 

In the following we prove that this conjecture is true under the 
hypothesis that the noise has a spherical distribution but the ra- 
dial density is not Lorentzian. More specifically, assuming only 
the stationarity of the noise process, it is proved that the con- 
densed density is asymptotically concentrated on the unit circle 
independently of the noise distribution. Moreover it is proved that 
if the noise distribution is spherical the condensed density Vn is 
independent of the specific spherical distribution, it is invariant 



oo 



f (^) = E a kZ 

k=0 



-k 
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by scaling, and, when n = 2, it is given by the uniform measure 
on the Riemann sphere. Furthermore in polar coordinates (/?, 9) 
the marginal condensed density w.r. to p (radial density) does 
not depend on 9 and the marginal condensed density w.r. to 9 
is uniform . Finally an integral representation of the condensed 
density is provided in the spherical case which shows that the 
radial density is not Lorentzian. 

The paper is organized in two sections. In the first one the prop- 
erties of the condensed density are studied. In the second one 
three numerical experiments confirming the derived properties 
are illustrated. 



1 The condensed density of the Pade poles 

Let us consider the transformation T = (T^\T^) that maps 
every realization a{p) of a = {a^, k = . . . , n — 1} to (£(0), 7(0)) 

given by a k {0) = E"4i TjMOM^j k = 0, . . . , n - 1, where G Q 
and Q is the space of events. It was proved in [31 Lem.2] that 
T(a) is defined and one-to-one a.e.. By noticing that the Cj are 
given by T^^(a), we define the condensed density as 



where the expectation is with respect to the density of a. 
Theorem 1 // {a^} is stationary 



Let us consider the periodic (circular) process {a^} obtained by 
repeating a finite segment of length p of the process {a^}. Then 



9 n/2 

h n (z) = -E E^-C) 



lim h n (z) = S(z — 1 




proof. 
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we have 

oo p—1 / oo 

f(z) = E *kz- k = E a J E 

fc=0 /i=0 \fc=o 

if |z| < 1 we get 

p-1 , z p 



-h 



f(z) = E a ^" „ , 

hence f (z) is a random rational function for every period p and 
its poles are the roots of unity. Therefore the Pade' denominator 
is a deterministic polynomial and the poles' condensed density is 
the counting measure on its roots. In the limit for p — > oo the 
periodic process {a^} converges in L 2 — mean to {a^} ([HI Sec. 
6]) and the counting measure on the roots of unity tends to the 
uniform measure on the unit circle. □ 

Theorem 2 If n = 2 and a = {9?[ao], S?[ao]? 3ft[ai], S[ai]} is 
spherically distributed with a density, then the poles' condensed 
density is the uniform measure on the Riemann sphere 

h 2 (z) 1 



7r(l + \z\ 2 ) 2 ' 

proof. 

If a n-dimensional random vector a is spherically distributed its 
joint density is ([8, th.2.9]) 

2 r (n/2) a 
/(Nl ) = ^rNI 9(\\a\\) 

where g(-) is the density of ||a|| and || • || denotes the euclidean 
norm. In the case considered n = 4. By making the change of 
variables T we get 7 = ao, C = an d 

lla|| 2 = |7l 2 (i + ICI 2 ). 

Noticing that the complex Jacobian of T is 7, the marginal on £ 
is then 

MO = hf(T(x))d-r = JL/ 7 g (|7lyTTKF)d 7 . 



as <p 
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By the change of variables 7 = 7^1 + |C| 2 ) and expressing the 
integral in real coordinates with real Jacobian I7I 2 we get 



hic) = *?l JTTW)W 9im TTW md ^ = 



' t^tP(|7|)^7^7- 



2^(1 +|C| 2 ) 2 



R 2 



But 



k I ^ff(l7l)«7<»7 = 1 



2 ^ 



because the marginals of a spherical density are spherical and we 
apply the formula above with n = 2 □ 

Theorem 3 If n > 2 and a = {3ft [a], S[a]} is 2n-variate spheri- 
cally distributed with a density, then the poles ' condensed density 
is the same independently of the specific distribution of a. More- 
over in polar coordinates (p, 9) the marginal condensed density 
w.r. to p does not depend on 9 and the marginal condensed density 
w.r. to 9 is uniform. Finally the condensed density is invariant 
by scaling a. 

proof. 

We first notice that the density of a is 



/(INI 2 ) = ^ll«l ll_2n 



27T f 



gum 



(2) 



where g(-) is the density of ||a|| because ||a|| = ||a||. We have 



K(z) = -E 



n 



n/2 

3=1 



=^r-i:/^-o)Nr- 2 "ff(ii«ii)<ifi 

n j=1<rn 
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and, by making the change of variables T, whose complex Jaco- 
bian is (0 Th.2]) 



^(c,7) = (- i r /2 n7in(o-a) 4 , 



n/2 



.7 = 1 j<h 



we get 



2 ^ r( 



n 



3-1 " 



n-l 

E 

fc=0 



n/2 

E 7/4 



2\ 1/2- 



5 



E 

fc=0 



n/2 

E 7/4 



^C^7 



/ y 



__ 2 ^ 2 r(n 
~ n ? tl ~2^ 







^(C (i U,7) 

.7-1 ^" (yn/2-l (pn/2 

^ < W2 2 \ 1/2 \ 

E 7aC* + 7i* 



/ n-l 

E 

^fc=0 



l,n/2 
E 7/4 + V 



2\ 1/2-n 
/ 



n-l 

E 

fc=0 



dC U) d 
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where = /i 7^ j} and 



J5(c uj ,^,7) = (-ir /2 n ih n (Cr-ar n (c-^r 



l,n/2 



/i=l r<h,rj^j 



Let us define 



Q7 = <2;(C (i) ,*) = XjXf e(E n ' 2 * n '\ X, G^ 2x ", Xj(h, k) = 41 



where 



x 



U) 
hk 



z k i, h = j 



(3) 



Then 



and 



n—l 

E 

fc=0 



l,n/2 

E 7*Cj + 7;*' 



i H Qa 



a j = i 



(4) 



where 



h f i\z)= J n (Cr-a) 4 n(Cr-^) 4 Ci(c (j) ^)^) 



and 



\h=l J 

We show now that Gj{C,^\ z) is independent of the specific form 
of the spherical density generator /(•). If 7 = {3ft [7], ^[7]} £ JR n 
and if Qj is the real isomorph of Qj, i.e. 



Qj = 



St[Qj], -$t[Qj] 



we have J_ H QjJ. = J^QiT- But then, by noticing that the real 
Jacobian of the transformation T is Jr = | Jc\ 2 , we get 

g^U) = / f'ff {?QnS ft - n 9 (if Qa) 1,2 )di 

M n \ h=1 J 

where Ah = I2 <S> e^e^ and is the h— column of the identity 
matrix of order n/2. Let us consider the transformation 

7 = QTi 



9 



which is well defined because Qj > and whose Jacobian is 
\Qj\ 1/2 . We then have 

r(n) r(n/2) , f 1 ^ 2 7 T Q7 1/2 ^Q7 V 1^1 f flr^ 1 -" „((*T*\V* 



r(n/2)W 2 |Q,|V 2 2W 2 ^ ^ ^ 
But then 



7 7 



7 T 7 9{{l T i) ) d l 



r(n/2) , /W 2 7 r QT 1/2 AQT 1/2 
2W 2 J M 



J 1 



7 7 



(y^Fi) #((i T i) 1/2 ) d i 



can be seen as the expectation of the function (n^i 2 



7 7 



of a n— dimensional spherically distributed r.v. 7. Therefore 

r(n) 



Gf.f/-(i) ^ = 

J - ' ; r(n/2)W 2 |Q J | 1 /2 

Let us define the matrix 



E 



l ff (l T Q- l,2 A h Q- ll2 l\ 

h=i 



1 T 1 



) \ 



$t[Qj }hh 

and prove that Bjh is idempotent. In fact 



Qj 1/2 A h Qj L/ ' z Qj L/2 A h Q] Lrz _ Q] l, ' z A h QfA h Q] i/2 _ 



-1/2 



-1/2 „ A-1/2 



1/2 



-1/2 



n>2 _ 



But 



0, e^ej 
7 



e/,ej, 
0, e^ej 
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because the diagonal elements of $s[Qj x ] are zero. Hence 



B 

We then have 



^•(C (i) , 



11 1 7 T 7 



(6) 



From [Bl Th.l.5.7,ii] it follows that w h = ( 7 ^ 7 ) has the beta 
distribution with parameters 1 and n/2 — 1 independently of the 
distribution of 7. As a consequence the distribution of h n (z) too 
does not depend on the distribution of 7. 

To prove the last part of the theorem, we notice that the spherical 
density (ED is invariant under the transformation 



The proof of Theorem 2 in [4] then holds and provides the re- 
quested results. □ 

Theorem 4 If n > 2 and & = {3ft [a], £y[a]} «s 2n-variate spher- 
ically distributed with a density, the poles condensed density is 
given by 



h n (z) = -\j- 2 J K Gn {z^ l) ) n ic.-ai 2 nic,-^ 

K Z7T ) 1 jff-2 r<h,(r,h)^l t+1 



where 



^-d9ft(C (1) )d3(C (1) * 



Eft 



s Pl (z,C 



2\ ™/ 2 + 



i^ G „(z,c (1) ) = ^ 

and y is a n-variate standard Gaussian random vector; 



l,n/2 

II (y r ^i h y 

h=l 
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sp 1 (z, C*- 1 )) are the Schur functions associated to the partition 

Pi = {ji, 32, ■ ■ Jn/2} 
which spans the minors of maximal order of X\ ( eq. \3^); 
sp^z,^ 1 ) are the Schur functions associated to the partition 

P2 = {jhh, ■ • - , in/2-1} 

which spans the minors of maximal order of the matrix obtained 
by Xi by canceling the h—th row. 



proof. 

From Theorem 3, without loss of generality, we can assume that 
7 has a zero-mean multivariate Gaussian density with identical 



covariance matrix. The ratio of quadratic forms 

can be rewritten as Wh = y T BjhY where y = jjX|. Therefore 
W/j, is a quadratic form in the variables y which are uniformly 
distributed on the n— dimensional sphere. But this distribution 
is spherical with characteristic function given by (P, 3.1.1] 

n 1 



*„(£) = oFig;Jl|l||j 



ttW 2 



In pZ3i Prop. 4] an explicit expression for Kc n (z, (^) = E 
when y is multivariate Gaussian is given, and it is claimed that 
this is also the result in the spherical case up to a constant which 
is a function of the characteristic function of the spherical distri- 
bution. It turns out that in the present case the constant is given 
by 

n n/2-l gr(,n/2+l)( ) T[n/2] 



(0) 



n/2 



2 n/2p[ 



n 



We then have 



E 



l,n/2 

n w ^ 

h=i 
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Moreover this expression is a symmetric function of z, C^. We 
then get by equations (Ej) and (EJ) 



\ Z7[ ) an/2-1 r<h,[r,h)^j r^j 



1 1,n / 2 

— n [Qj'Udc^ . 



\Qj\ 1/2 



h=l 



Remembering that the real Jacobian of the transformation T is 
Jr = \Jc\ 2 , we have 



, 9 ^72 / ^„(^c (i) ) n Kr - ai 8 n ic. 

l Z7r J ron-2 r<h,(r,h)^j r^j 



\Qj 



h=l 

because \Qj\ = \Qj\ 2 . But, by Binet-Cauchy formula, we have 

2 



jl<?2<-<jn/2 



[l 2 . 


..n/2^ 




• • Jn/2 J 



where Xn 



( \ 2 ...n/2^ 



is a minor of maximal order n/2 of Xj 



J2 • • • y 

From [12] we have 



1 2 ...n/2 

^jl J*2 ■ ■ • jn/2 J 



s Pl (zx ij) ) n (Cr-a)n(Cr-^) 

r<h,(r,h)^j r^j 



where sp^z, Q^) is the Schur function associated to the partition 
Pi = {ji, j2, . . . , Jn/2}? which is a symmetric polynomial with 
positive integer coefficients (Jack function with a = 1). Hence 



1 = n ICr 


- Ch\ 2 n Kr - 


;| 2 £ 




r<h,(r,h)^j 




h<32<-<3n/2 
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But [Q7V = [ad jgi ?)]/ " and 



[adj(<9j)]fch = E 

3l<h<-<3n/2-l 



1, 2, .../i-l,/i + l...n/2 ^ 



Jn/2-1 y 



p 2 



n Kr-ai 2 n 

r<k,(r,k)^j,h r^j^h 



where P 2 = {ji, J2, • • ■ , in/2-1} and 



hg>(z) = 



(2tt 



1 Rn-2 r<h,(r,h)^j r^j 

l,n/2 f 

n n ic-ai 2 n iCr-* 

fc=l \r<k,(r,k)^j,h r^j^h 

2 



6— n 



n l,n/2 / v 
llft=l Lp 2 



sp 



•^,C (1) ) 



2\ «/2- 



f^(C (i) )^(C (i) )- 



But it turns out that 

l,n/2 / 

n n ic.-ai 2 n iCr 

/i=l \r<k,(r,k)^j,h r ¥ z j¥ z h 

therefore 



n ICr— ar -4 n ic-^ 

r<h,(r,h)^j r^j 



n— 4 



(2tt)"/ 2 



ntf Ep 2 



(ep x 



n ic,-ai 2 n ic,-^i 2 - 

r<h,(r,h)^j r=fj 

sp 2 (z,& h) ) 2 ) „ 

^(C (J) )^(C (J) )- 



sp 



i(-,C (1) ) 



2\ n/2+1 



Because of the symmetry of the Schur polynomials and of K n (z, 
all h%\z) : j = 1, . . . , n/2 are equal. Therefore by equation (Hj) 
we get the thesis. □ 
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From the result of the theorem it follows that, expressing h n (z) 
in polar coordinates and taking the marginal with respect to the 
modulus (radial density), we get a much more complicated ex- 
pression than a Lorentzian function, thus disproving a part of the 
conjecture mentioned in the Introduction. We also notice that 
evaluating h n (z) is not an easy task because computing Schur 
functions is far from trivial (0). 

2 Numerical examples 

In order to illustrate the results obtained in Theorem 3 and 4 
three numerical experiments were performed. In the first one a 
sample of cardinality 2 • 10 6 of a multivariate complex a— stable, 
centered, symmetric colored noise with a = 0.5 and n = 4 was 
generated. The Pade' poles were computed as well as the empiri- 
cal density of their modulus. In fig. 1 this estimate is represented 
by dots. Then the integral representation of the poles condensed 
density derived in Theorem 4 was transformed in polar coordi- 
nates z = pcos(6) with Jacobian p. After Theorem 3 we know 
that the phase distribution is uniform. Therefore the radial dis- 
tribution is obtained by multiplying by 2irp the integral repre- 
sentation given in Theorem 4. The integral was approximated by 
numerical quadrature where, in this case 

K Gi (z, = tr(5i)tr(5 2 ) + 2tr(5 1 £? 2 ), 

and is plotted in fig.l as a solid line. We notice that the fit is 
quite accurate. 

In the second experiment a 4-variate complex Gaussian white 
noise was generated. The same computations as above were per- 
formed and the results are plotted in fig. 2 confirming that the 
condensed density in the Gaussian case is the same than in the 
a— stable one as claimed in Theorem 3. 

In the third experiment a sample of cardinality 2 • 10 6 of a 6- 
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variate complex Gaussian white noise was generated. The same 
computations as above were performed. In this case 



K G6 (z, C W ) = tr(Bi)tr(B 2 )tr(B 3 ) + 2 [tr(5i5 2 )tr(5 3 ) + 
tr(5 1J B 3 )tr( J B 2 ) + tr(5 2 S 3 )tr(5i)] + 8tr( J B lJ B 2 5 3 ). 

The results are plotted in fig. 3. Also in this case the fit is quite 
accurate. 
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Fig. 1. The poles condensed density for n = 4, computed by numerical integration 
(Th.4) (solid) and the empirical density based on 2 • 10 6 replications of a colored 
noise with a symmetric, centered a— stable density with a = 0.5 (dotted). 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Fig. 2. The poles condensed density for n = 4, computed by numerical integration 
(Th.4) (solid) and the empirical density based on 2 • 10 6 replications of a complex 
Gaussian white noise (dotted). 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Fig. 3. The poles condensed density for n = 6, computed by numerical integration 
(Th.4) (solid) and the empirical density based on 2 • 10 6 replications of a complex 



Gaussian white noise (dotted). 
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